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I. INTRODUCTION 


Elastic and inelastic collisions between atoms and molecules and/or between molecules 
and molecules are of great scientific interest in the fields of physical-chemistry and chemical- 
physics. The reason is that such processes can provide valuable information about fundamen¬ 
tal interactions between different chemical species, their chemical properties, their energy 
transfer quantum dynamics and many other properties. The pioneering studies focussed on 
light atoms and molecules, because of their simple nature. For three- and four-atomic sys¬ 
tems with a small number of electrons, the potential energy surface (PES) can be computed 
with relatively high accuracy^”-. Consequently, for these systems one can then more easily 
test different dynamical methods, such as, classical, semi-classical, quasi-classical trajectory, 
and pure quantum-mechanical computational formulations and compare the results with 
available experimental data in a controlled fashion. The test methods devised could then be 
applied to more complex many-atomic systems, wherein a controlled comparison is not pos¬ 
sible. Among these small systems the four-atomic H2-I-H2 and H2-I-HD scattering processes 
have attracted signihcant attention not only in chemical physics but also in astrophysics. In 
astrophysics H2 and HD play an important role, because of their abundance in the molec¬ 
ular cloud of the universe^*^. Together with the H2 -|- H2 collision, the HD -|- H2 collision 
is also of signihcant importance in the astrophysics of the early universe. Specihcally it is 
important in the modeling of pre-galactic clouds and planetary atmospheres; in the cooling 
of primordial gas and in the formation of stars^.^— . 

In Ref.— a semiclassical treatment of H2-H2 scattering is formulated. In Ref.— the author 
developed and applied a rigid rotor model to study rotational excitation in H2-H2 by applying 
quantum close-coupling scattering calculations. In this approach the distance between the 
hydrogen atoms in both H2 molecules was hxed at a constant value based on avarage. This 
model was applied to many different atomic and molecular systems, see, for examplei^'i^'"— . 
The main goal of the work— was to compute rotational thermal rate coefficients in the 
H2-I-H2 system at low temperatures of astrophysical interest. Quantum-mechanical close- 
coupling calculations for three-dimensional collisions of para-H2 and ortho-Y{2 with HD are 
performed in Ref.— where the HD-H2 potential is derived from the H2-H2 potential. 
A quantum dynamical study of H2-H2 collisions is reported for both ortho- and para-B.2 in 
Refs.— ^— . In Ref.—, the authors considered H2-I-H2 and D2-I-D2 rotational inelastic scat- 
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terings with the use of the H2-H2 potential energy surface (PES) from Ref.-. In Ref.— a 
full six-dimensional scattering computation has been performed taking into account vibra¬ 
tional relaxation in the H2-I-H2 collision. In that study the H2-H2 PES from Ref.- was used 
and its anisotropy properties have been studied at low temperatures: 20 K < T < 300 K. 
HD-I-HD scattering has been studied experimentally in Ref.— and theoretically in 2 ^ for a 
wide range of collision energies. A comprehensive computational and experimental study of 
total cross section in H2-H2, D2-D2, and HD-HD scattering for both ortho and para H2 and 
D2 has been reported in Ref.—. Measurements of energy transfer rates in HD-I-HD^S and 
H2+H2^>^ collisions have also been performed. 

However, realistic theoretical investigations of the low-energy HD-I-H2 collision are lack¬ 
ing, although preliminary quantum calculation of this process has been reported in Refs.— 
Schaefer— calculated rate coefficients for the excitation of HD by H2, for the low-lying rota¬ 
tional levels using a modihed older potential for HD and H2. Flower— performed an improved 
calculation of HD-H2 scattering with Schwenke’s H2-H2 PES using a larger rotational basis 
set^. In this paper we report an improved calculation of this problem using a realistic HD-H2 
PES derived from Hinde’s recent H2-H2 PES^. In two papersi^^l^ the PES from work- has 
been applied together with a pure quantum-mechanical dynamical approach. The surpris¬ 
ing thing is, that the results of worksl^ii^ are closer to the results of older work— than to^. 
Therefore, there is a need to carry out new computations with newer PES between HD and 
H2. 

While the exchange symmetry is broken in HD-I-H2 it still posses many similarities with 
the H2-I-H2 system. The PESs of H2-H2 and HD-H2 should basically be the same six¬ 
dimensional function. This fact follows from the general Born-Oppenheimer approach^. At 
the same time the two collisions: H2-I-H2 and H2-I-HD, should have rather different scattering 
outputs. This is because the H2 and HD molecules have fairly different rotational constants, 
internal symmetries and as a result a different rotational-vibrational spectrum. The HD-H2 
PES can be derived from the H2-H2 PES by shifting the center of mass (c.m.) of the H2 
molecule to the c.m. of the HD molecule. Once the exchange symmetry is broken in H2-H2 by 
replacing the H with the D atom in one of the H2’s then one has the new HD-H2 PES. In this 
fashion, we constructed the HD-H2 PES from the H2-H2 PES of Hinde^ employing all parts 
of the full HD-H2 interaction including the HD’s dipole moment. Using this HD-H2 PES 
we carried out pure quantum-mechanical calculations for inelastic collisions of rotationally 
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excited HD and H2 molecules, i.e. the process: 

HD(jO + H2(j2) -> HD(jO + R.if,). ( 1 ) 

The scattering cross sections and their corresponding thermal rate coefficients are computed 
using a non-reactive quantum-mechanical close-coupling approach. The four-atomic system 
is shown in Fig. [H 

In Sec. II we briefly describe the quantum-mechanical approach used in this paper. Sec. 
Ill includes the computational results. We compare the cross-sections and rates with those of 
other authorsiiii^, and our previous calculationsi^ii^, where a different HD-H2 PES, derived 
from the well-known Boothroyd-Martin-Keogh-Peterson (BMKP) H2-H2 PES^, was used. 
Discussion and conclusions are provided in Sec. IV. The corresponding procedure to obtain 
a modihed HD-H2 PES from the existing H2-H2 surface^ is presented in Sec. |Vl Atomic 
units {e = nie = h = 1 ) are used throughout this paper. 


II. QUANTUM-MECHANICAL APPROACH 


In this section we provide a brief account of the present quantum-mechanical close¬ 
coupling approach following the method in Ref.—. The HD and H2 molecules are treated 
as linear rigid rotors. The model has been applied in few previous works^^ii^. In all our 
calculations with this potential the bond length was hxed at 1.449 a.u. or 0.7668 A for the 
H2 molecule and 1.442 a.u. for HD which is 0.7631 A. The Schrodinger equation for the 
( 12 ) -|- ( 34 ) collision in the c.m. frame, where ( 12 ) and ( 34 ) are diatomic molecules formed 
by atoms 1-4 is^^: 


r Pi 
_^ 

2 M 
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^R 2 


+ V (Ri, i?2, R 3 ) — E 


12 


T(Ai, A2, R3) — 0 . 


( 2 ) 


2 jji\R^ 2H2R2 

Here is the momentum operator of the kinetic energy of the collision, P3 is the collision 
coordinate, whereas Pi and P2 are relative vectors between atoms in the two diatomic 
molecules as shown in Fig. [H and are the quantum-mechanical rotation operators of 

the rigid rotors, fii and ^2 are the reduced masses of the HD and H2 molecules and AI12 is the 
reduced mass of the two molecules. The vectors Pi(2) are the angles of orientation for rotors 
(12) and ( 34 ), respectively; V{Ri, R2, R3) is the PES of the four-atomic system ( 1234 ), and 
E is the total energy in the c.m. system. The use and modihcation of the original H2-H2 
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PESs V{Ri, R2, R3) is discussed in Appendix A, i.e. Sec. 5 . The cross sections for rotational 
excitation and relaxation can be obtained from the A-matrix S^^,. The cross sections for 
excitation from HD(ji,mi)+H2(j25’^2) to HD+112(^2^2) summed over the hnal 
angular momentum projections {m'^m'2) and averaged over the initial projections (mim2) of 
the HD and H2 molecules of angular momenta ji and j2 are given by: 




TT 

(2jl + l)(2j2 + 1)^00' 


(2J + 1)|(5oq' 


siAE)(^ 


( 3 ) 


The kinetic energy is e = E — + 1 ) — B2j2{j2 + 1 ), where Bi = 44.7 cm“^ 

and B2 = 60.8 cm“^ are the rotation constants of rigid rotors ( 12 ) and ( 34 ) respec¬ 
tively, they are shown in Fig. [T] Next, E is the total energy of the system, J is total 
angular momenta of the four-atomic system, a = (jij2ji2T), where ji -|- j2 = ji2 and 
ji2 + L = J, kaa' = ^2^^12(77 + Ea — Ea>) is the channel wavenumber and Ea(a') are 
rotational channel energies. Finally, the relationship between the rotational thermal-rate 
coefficient at temperature T and the corresponding cross section 

can be obtained through the following weighted average formula: 
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( 4 ) 


where is Boltzman constant and e = E — Ej^ — Ej^ is pre-collisional translational energy 
at the translational temperature T, ks is Boltzmann constant and Sg is the minimum kinetic 
energy for the levels ji and j2 to become accessible. 


III. NUMERICAL RESULTS 

We used the MOLSCAT progrant^i to solve the Schrodinger Eq.([ 2 ]). Convergence was 
obtained for the integral cross sections, cr(jj, Jii2,£), with respect to the variation of the 
variables utilized in all considered collisions at different collision energies. For the inter- 
molecular distance R3 we used from R^min = 3.0 a.u. to R^max = 30.0 a.u. We also applied 
a few different propagators included in the MOLSCAT computer program, and our calcu¬ 
lation show that D. Manolopoulos’s hybrid modihed log-derivative propagator technique^ 
would be quite numerically stable and a time effective approach. This method is used in the 
majority of the calculations. 
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The maximum value of the total angular momentum J was set at 44 while the number of 
levels Nivi in the basis set of HD + H 2 was set at 42. Specifically, in the case of HD(ji)+para- 
H 2 (j 2 ), ji it has values 0, 1, 2, and 3 and 72 has values 0, 2, and 4. This combination generates 
the total number of included levels iV;^/=34. In the case of HD+orfho-H 2 the parameter ji 
has values 0, 1, 2, and 3 and 72 has values 1, 3, and 5. This results in the total number 
of levels Nm=A2 . A number of test computations with higher values for the ji and j 2 
parameters have been carried out. For example, in the case of HD+para-H 2 , ji was taken 
as 0, 1, 2, 3, 4 and 72 as 0, 2, 4, 6. This ji/j 2 combination produced Nm = 74. We obtained 
similar results in both cases, confirming the convergence of the calculations. 

Because the HD+H 2 total rotational energy transfer cross sections Q has shape reso¬ 
nances at low energies a large number of energy points are needed in order to effectively 
reproduce them. We used up to 250 energy points in each computation for each specihc ro¬ 
tational transition in the HD and H 2 molecules considered. More space discretization points 
were used at low collision energies and fewer points in the higher-energy sector. 

Below we present results for cross sections and thermal rate coefficients for different 
quantum-state transitions in HD and H 2 molecules. We reproduced shape resonances in the 
low velocity region, which are very important in the cooling of the astrophysical media at low 
temperatures. We compared them with the older quantum-dynamical results of Schaefer— 
and Flower— . We also present resultsl^i^ using a PES obtained from a modification of the 
BMKP H 2 -H 2 PES^, viz. Sec. |V] (Appendix A). In Figs. [2]|3] and l4l we show four different 
results for the integral cross sections in the collisions: 

HD(1) + H2(0) ^HD(0) + H2(0), (5) 

HD(2) + H2(0) ^HD(1) + H2(0), (6) 

HD(2) + H2(0) ^HD(0) + H2(0), (7) 

respectively. One can see that the new cross sections obtained with the present PES have the 
same structure and shape, but also have substantially larger (~ 60%) values at medium en¬ 
ergies {v > 200 m/sec) when compared with the results obtained from the modihed BMKP 
PES^ and the older result from Schaefer’s calculations!^. Processes dH]) and m are also 
important collisions from the astrophysical point of view. They represent transitions from 
the ji = 2 and j 2 = 0 state of the HD(ji)-|-H 2 (j 2 ) system. It is seen from Table [T] that 
the corresponding rotational energy of the system is 268.2 cm“^ and rotational relaxation 
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processes to states of lower energy should be relevant. The cross sections from Ref.— are 
available only for higher collision velocities, i.e. v > 300 m/s. The present cross sections 
have the same behavior as previous results, but they have larger values than the old ones, 
specially at higher energies. Next in Fig. |5]we exhibit the thermal rate coefficients corre¬ 
sponding to the processes considered above. Here, we additionally include the older results 
from Ref.—, where Schwenke’s H 2 -H 2 PES^ was used. One can see that the new rates are 
substantially larger than those of other calculations. In general, it means that the contribu¬ 
tion of the HD-I-H 2 collision to the HD-cooling function can have even larger contributions 
than previously expected. 

Now we consider the inelastic cross sections in HD-I-H 2 for higher rotational energies. For 
the systems HD(0)-|-H2(2) and HD(1)-|-H2(2), one can see from Table [I] that the correspond¬ 
ing rotational energies are 364.8 cm“^ and 454.2 cm“^. In Fig. |6] (upper plot) we show the 
integral cross sections for the process: 

HD(0) + H2(2) ^HD(2) + H2(0). (8) 

While we obtain a relatively good agreement between cross sections calculated with the older 
modihed BMKP PES and Hinde’s PES, there is a dramatic difference with the corresponding 
result from Ref.—. The present cross sections are larger than those of Ref.— by a few orders 
of magnitude. In Fig. IH] (lower plot) we show the results for the corresponding thermal rate 
coefficient k 2 o^o 2 (T). It is seen, that the present rates and those obtained by Flower— have 
a flat temperature dependence, whereas the rate of Schaefer—, although smaller than other 
results, increases monotonically with energy. Because the thermal rate of reaction ([8]) is 
relatively large, one can conclude that this channel can also make a substantial contribution 
to the astrophysical HD-cooling function. 

In Figs. ITlfTOl we show results for the total cross section in the inelastic scattering from 
the state HD(1) -I- H2(2) with rotational energy 454.2 cm“^, viz. Table HI All de-excitation 
processes have been computed for this state, namely: 


HD(1) + H2(2) HD(0) + H2(2), 

(9) 

HD(1) + H2(2) HD(2) + H2(0), 

(10) 

HD(1) + H2(2) HD(1) + H2(0), 

(11) 

HD(1) + H2(2) ^ HD(0) + H2(0). 

(12) 
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The present cross sections exibit fairly good agreement with the results computed with the 
modified BMKP PES^, and also with older results from Ref.—. The general shape and 
trend of the behavior of these cross sections are the same in all cases. Also, there is a 
relatively small bump in the cross sections of the processes f[TU]) - f[T^ at collision velocity 
~1100 m/sec which is also reproduced by two PESs. One can see that the process (E]) can 
make a signihcant contribution to the total HD-cooling function because its cross section is 
rather large relative to other channels. Finally, in order to carry out new computations of 
the astrophysical cooling function. Table Hi] includes the relevant thermal rate coefficients 
for the HD + para-H 2 case in the temperature range from 2 K to 300 K. Next, as a test, 
we choose the initial state HD(2) + H2(2) with a relatively higher total rotational energy: 
633 cm“^. The integral de-excitation cross sections from this state to different lower energy 
rotational states are shown in Fig. [TT] One can see a fairly good agreement in the shape 
of the curves between various rotational transition results for the integral cross sections 
computed with the two different PESs. 

In Fig. [12] we show three different rotational transition cross-sections for the ortho¬ 
hydrogen case. Here we chose low lying rotational levels of the two molecules: HD(1) 
-I- H2(1) and HD(2) -|- H 2 (l). Some results from older worksii'i^'I^ are also presented in the 
figure together with results computed with the newer modihed PES from^. The correspond¬ 
ing thermal rate coefficients are shown in Fig. [13] where the older results from Refs.— 
are also presented for comparison purposes. Figures [T4| and [I5] include results for thermal 
rate coefficients for the transitions from higher rotational states, e.g., HD(0) -|- H2(3) and 
HD(1) + H2(3). 

Finally, in Sec. I VI I (Appendix B) we present our new results for the thermal rate co¬ 
efficients which can be used in subsequent computation of the astrophysical HD-cooling 
function. In Table m we show thermal rates of different de-excitation processes in low- 
temperature HD -|- para-B .2 rotational energy transfer collisions and in Table IIIII the same 
data for the HD -I- ortho-B .2 case. 


IV. SUMMARY AND CONCLUSIONS 


State-to-state close-coupling quantum-mechanical calculations for rotational de-excitation 
cross-sections and corresponding thermal rate coefficients of the HD-|-o-/p-H 2 collisions are 




presented using a linear rigid rotor model for the HD and H 2 molecules. The symmetrical 
H 2 -H 2 PES of Ref.- has been appropriately adopted for the current non-symmetrical HD+H 2 
system by appropriate translation and rotation. These geometrical operations lead to a new 
set of angle variables 6 [, and <^2 fhe Jacobi few-body coordinates, a new length of the 
intermolecular distance R 3 and, as a result, to a new HD-H 2 PES. For comparison purposes 
in this paper we carried out a few calculations with the use of the older BMKP PES^ which 
was also modihed for HD-H 2 . A test of convergence and the results for the cross-sections 
with the two PESs are obtained for a wide range of values of different parameters. It is seen 
from Figs. [2113] and 01 that in most cases for the lower number quantum transition states the 
rotational energy transfer cross sections obtained with the use of the modihed Hinde’s PES 
have higher values than the cross sections computed with the use of the modihed BMKP 
PESl^J^. The same situation occurs in comparisons within the older work—. In Fig. [5] the 
corresponding thermal rate coefficients are presented. In this case we include results from 
the paper— too. As can be seen the new thermal rates (solid lines) have substantially higher 
values than the other results. For example, in the rotational energy transfer, process ([5]), 
new thermal rates are ~2 times higher than the other corresponding rates. Because ([5]) is 
considered as one of the main contributors to the resulting HD cooling process, one can say 
that the new HD-cooling function may have substantially greater values than the previous 
calculations!^. Next, Fig. [ 6 ]represents results for a very interesting process (| 8 ]). The interest 
in this channel lies in the fact that both molecules change their rotational quantum numbers 
by Aj=2, therewith HD becomes excited and H 2 de-excited. In this case we obtained a 
significant deviation from the results of work—. However, the cross sections obtained with 
the modihed versions of the Hinde and BMKP PESs are fairly close to each other and have 
fairly large values. Therefore, this process probably can make a contribution to the HD 
cooling process. In addition we would like to note, that in work— the process has also been 
computed and discussed, in which signihcant differences from!^ were also found. 

Further, the following four Figs. (j7 |) - fiT0|) represent our integral cross sections and corre¬ 
sponding results from work— for the following de-excitation collisions: ([2]), fllUp . fITT]) . and 
(1T2|) . In all of these processes the initial state of HD has the rotational state ji=l and H 2 
has the rotational state j 2 = 2 . The corresponding total rotational energy of the molecules 
can be found in Table [H The results in Fig. ([^ are in fairly good agreement with each other 
and have relatively large values, therefore this specific channel could also make a substantial 
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contribution to the cooling function. It was found, that in Figs. fl 8 |)- (fT 0 D both surfaces, 
namely, Hinde’s and BMKP both accurately reproduce the position of a shape resonance at 
collision velocity v ~1300 m/sec. Fig. ITT] shows the resulting de-excitation cross sections 
from the highly located rotational quantum level. Finally, our analysis in this paper would 
not be complete if we did not undertake computations for the ortho hydrogen case as well. 
Fig. [12] shows results for the lower lying rotational quantum numbers of HD, namely, ji=0 
and 1, and H 2 : j 2 =F It is seen that in this case the cross sections obtained with the Hinde 
PES are signihcantly larger than other results. The corresponding thermal rate coefhcients 
are presented in Fig. [T31 Again, as in previous par a-hydrogen cases the new rates obtained 
with Hinde’s potential are larger than previous results. At the low density limit and taking 
into account the critical density concept the total cooling function can be computed with 
the use of the following formula: 


Ahd(T) = '^hd(ji)'^H2 {j^)^jij2^j[j2 (13) 

which is in units of [ergxcm“^xs“^]. Here, is the emitted photon energy, 

it the thermal rate coefficient (| 1 |) corresponding to the rotational transitions 
jij 2 Therefore, increasing the knowledge of rotational and possibly vibrational 

excitation and de-excitation rate constants, in atomic and molecular hydrogen- 

hydrogen collisions, such as HD/H 2 -I-H 2 , HD/H 2 -I-H etc, is important in order to under¬ 
stand and be able to model the energy balance in the interstellar medium. For comparison 
purposes it would be very useful and interesting to carry out new computations of the 
rotational-vibrational integral cross sections and corresponding thermal rate coefficients for 
a low-energy HD-I-H collision. In this case a different H 3 PESs from papers^*^ could be 
applied. 
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V. APPENDIX A: HD-H2 POTENTIAL ENERGY SURFACES 


A few important modifications to the Hinde H 2 -H 2 PES^ were needed for the current 
non-symmetrical four-atomic collision ([T]). The application and modihcation of the original 
H 2 -H 2 BMKP PES were published in Ref.—. Below in this paragraph we briefly describe 
the procedure. To compute the distances between the four atoms the BMKP PES uses 
Cartesian coordinates. Consequently, it was necessary to convert spherical coordinates used 
in the close-coupling method to the corresponding Cartesian coordinates and compute the 
distances between the four atoms followed by calculation of the PES^^ii^. This procedure 
used a specihcally oriented coordinate system OXYZ. As a hrst step we needed to introduce 
the Jacobi coordinates {Ri, R2, R3} and the radius-vectors of all four atoms in the space- 
hxed coordinate system OXYZ: {fq, P 2 , fh, u}- Then the center of mass of the HD molecule 
has been relocated at the origin of the coordinate system OXYZ, and the R3 was directed to 
center of mass of the H 2 molecule along the OZ axis. Thus, one could obtain the following 
coordinate relationships: R 3 = {i? 3,03 = 0 , $3 = 0 }, with ©3 and $3 the polar and 
azimuthal angles, Ri = fq — fq, R 2 = fq — fq, fq = ^Ri and rq = (1 — where 

^ = m 2 /{mi + 772 , 2 )^. Further, we adopted the OXYZ system in such a way, that the 
HD inter-atomic vector Ri lies on the XOZ plane. Then the angle variables of Ri and R2 
are: Ri = {0i,$i = vr} and R 2 = {© 2 ,^ 2 } respectively. One can see, that the Cartesian 
coordinates of the atoms of the HD molecule are^: ?q = {xi = sin ©i, 7/1 = 0, Zi = 
.^RiCosOi}, fq = {x 2 = -(1 - sin © 1 , 1/2 = 0,^2 = -(1 -O-^iCosOi}. Dehning ( = 
{ 1713 + 1714 ), we have fq = R 3 — {1 — C)R 2 , fq = R 3 +CR 2 , and the corresponding Cartesian 
coordinates are: fq = {xs = —(1 — C)i ?2 sin ©2 cos $ 2 , Vs = —(1 — 0^2 sin ©2 sin $ 2 , Z 3 = 
R 3 - {1 - C)R2Cos 02}, fq = {xi = C-R 2 sin© 2 Cos$ 2 , Ua = C-R2sin02sin<h2, Z 4 = 
R 3 + CR 2 cos © 2 }. In such a manner the cartesian and the Jacobi coordinates are represented 
together for the four-atomic system HDj-H^I^. 

The Hinde H 2 -H 2 PES^ is a six-dimensional surface which was constructed using recent 
Raman spectrum data of the (H 2)2 dimer and which accurately describes the dimer’s van 
der Waals well^. It was demonstrated that this PES gives IR and Raman transition energies 
for the {para-B. 2 ) 2 , {ortho-D 2 ) 2 , and (para-H 2 ) —(or/:ho-D 2 ) dimers and is in good agreement 
with experimental data. 

The method to make the Hinde H 2 -H 2 PES suitable for the non-symmetric system HD-I-H 2 
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is based on a geometrical transformation technique, i.e. a rotation of the three-dimensional 
space and the corresponding space-hxed coordinate system OXYZ. The new global PES 
depends on six variables (Jacobi coordinates) — |-Ri|, |^ 2 |, I-R 3 I, 6 ^ 1 , 6 * 2 , and <^2 — as shown 
in Fig. [T] together with corresponding quantum angular momenta. The initial geometry of 
the system is designed in such a way that the Jacobi vector connects the c.m.’s of the two 
H 2 molecules and is directed over the OZ axis. We laid out OXYZ in such a manner that 
the Jacobi vector Ri lies in the XOZ plane. The vector R 2 can then be directed anywhere. 
The spherical coordinates of the Jacobi vectors are: Ri = R 2 = {i? 2 , 6 ^ 2 , 932 }, 

and R 3 = {i? 3 , 0, 0}. Because we used the rigid rotor model, the lengths of the HD and H 2 
molecules are hxed at equilibrium values, e.g., i?i=0.7631 a.u. and i?2=0.7668 a.u., thus 
leaving us with four free variables. We replace one hydrogen atom H with a deuterium atom 
D, thus shifting the c.m. of one H 2 molecule to another point, that is from O to O'. The 
length of the vector x is x = |i?i|/ 6 . This is seen in Fig. [161 Then we shift the original 
coordinate system OXYZ along the vector Ri to the new one O'X'Y'Z'. The origin of the 
new system, i.e. O', lies on the c.m. of HD. A new Jacobi vector R!-^ is defined connecting 
the c.m.’s of the HD and H 2 molecules. The new intermolecular distance between HD and 
H 2 is: 

= \JR\ + x'^ — 2 xi ?3 cos 6 * 1 . (14) 


Now, if we rotate O'X'Y'Z' around its O'Y' axis in such a way that the OZ' axis is directed 
over the vector we obtain a new coordinate system O'X"Y"Z" which should be well 
designed to carry out computations for the HD-I-H 2 collision. The rotational angle rj satisfies: 


_ ^2 


COST] = 

sinr; = 




3-^3 


X 


R'n sin 9i 


(15) 

(16) 


This transformation converts the initial Jacobi vectors in OXYZ: Ri = {i?i, 6 ^ 1 , 0 }, R 2 = 
{R 2 , O 2 , 9 ? 2 } and R 3 = {i? 3 , 0 , 0 } to the corresponding Jacobi vectors with new coordinates 
in the new O'X'Y'Z': R'^ = {R'^, (p' 2 }, R'^ = {R' 2 , 9'^, 0} and R'^ = 0, 0}. 

The coordinate transformations from OXYZ to O"X"Y'Z" changes the original Hinde’s 
PES to the new HD-H 2 PES. Rotation of the coordinate system from O'X'Y'Z' to 0'X"Y"Z" 
results in a corresponding transformation of the coordinates of the 4-body system as well 
as changing the distance between the two molecules. One then has the following relations 
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between new and old variables^: 


d'^ = arccos(cos Q\ cos rj + sin 9i sin rj ), 


0 ' 

^2 


/ cot6>isin?7\ 

: arccos cot 02 cos rj H-^- 

V sin 02 / 


( 17 ) 

( 18 ) 
( 19 ) 


In the calculation of HD+H2 with Hinde’s PES one has to use new coordinates 9[, 6*2, ^2, -R3. 
However, the original potential has been expressed through the initial H2-H2 variables, i.e. 
9i, 92, $2 and R^. Hence they have to be transformed using flT7)) - f[T^ . Therefore, in the 
case of the non-symmetrical HD+H2 collision one should use the formulas ffT 7 )) - ffT^ together 
with (lT^ - (fT6|) and the expression (1T4)) for the new distance R'^ between the center of masses 
of the H2 and HD molecules. 

In general, any consideration of the HD+H2, D2+HD or D2+D2 systems should begin with 
the original H2-H2 PES. This six-dimensional function comprises a symmetrical surface over 
the OZ coordinate axis. This is shown in Figs. 1 and 16 . In spherical coordinates the surface 
can be described by six variables: Ri, 9i, R 2 , 6*2, R 3 , and ip 2 . The variables are also shown 
in the figures. The H2-H2 PES was obtained in the framework of the Born-Oppenheimer 
model^ and can be considered as a symmetrical interaction field. When considering non¬ 
reactive scattering problems with participation of hydrogen molecules one needs to solve the 
Schrodinger Eq. (| 2 ]) with the H2-H2 potential V{Ri, R 2 , R 3 ). The solution/propagation runs 
over the R3 Jacobi vector (please see Fig. 16 ). Therefore, in the case of the symmetrical 
H2-I-H2 and D2-I-D2 collisions one can use the original H2-H2 PES as it is, i.e. without 
transformat ions. 

However, in the case of the non-symmetrical (or symmetry-broken) HD-I-H 2 /D 2 or 
HD-I-HD scattering systems one should also apply the original H2-H2 interaction field 
(PES), but the propagation (solution) of the Schrodinger equation runs, in this case, over 
the corrected Jacobi vector R'^ which is directed over the new OZ" axis, as is shown in Fig. 
16 . 


VI. APPENDIX B: HD+0-/P-H2 THERMAL RATE COEFFICIENTS 

New data for the thermal rate coefficients, Eq. (jl]), are listed below. The 

results are obtained with the modified Hinde H2-H2 PES^: Table UTI includes (T) at 
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low temperatures for the HD + para-H 2 rotational energy transfer collisions and Table IIIII 
the same information for the HD + ortho-Y \.2 case. The astrophysical HD-cooling function 
can be computed with the use of formula (IT^ . 
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TABLE I. Total rotational energy £rot in the four-atomic system: HD(ji)+p-H 2 (j 2 ) and HD(ji)+o- 
H 2 (j 2 )- Here, £rot = -|- 1) + B 2 j 2 {j 2 + 1)) where Bn 2 ) are the rotational constants of rigid 

rotors ab and cd respectively. 


£rot (cm“^) 

HD(ji) 

P-B2{j2) 

£rot (cm“^) 

HD(ji) 

0-B.2{j2) 

0.0 

0 

0 

121.6 

0 

1 

89.4 

1 

0 

211.0 

1 

1 

268.2 

2 

0 

389.8 

2 

1 

364.8 

0 

2 

658.0 

3 

1 

454.2 

1 

2 

729.6 

0 

3 

536.4 

3 

0 

819.0 

1 

3 

633.0 

2 

2 

997.8 

2 

3 
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TABLE II. Low temperature rotational de-excitation thermal rate coefficients kij^iiji{T) in the 
HD(i) +para-H 2 (j) —>■ HD(z') -|-para-H 2 (/) collision. All results are multiplied by a constant value 
a = 10^^. The data are in the unit cm^ sec“^. 


r(K) 

10-00 

20-10 

20-00 

02-20 

02-10 

02-00 

12-02 

12-20 

12-10 

12-00 

2 

4.774 

6.892 

1.115 

0.914 

3.88E-02 

2.45E-03 

4.960 

2.22E-02 

5.61E-03 

6.01E-04 

4 

4.099 

4.817 

0.670 

0.553 

2.31E-02 

1.46E-03 

4.194 

1.74E-02 

4.40E-03 

4.56E-04 

6 

3.780 

4.007 

0.485 

0.444 

1.83E-02 

1.16E-03 

3.841 

1.55E-02 

3.96E-03 

4.04E-04 

8 

3.574 

3.580 

0.390 

0.393 

1.60E-02 

l.OlE-03 

3.617 

1.44E-02 

3.74E-03 

3.76E-04 

10 

3.431 

3.321 

0.335 

0.363 

1.47E-02 

9.26E-04 

3.464 

1.38E-02 

3.61E-03 

3.60E-04 

12 

3.329 

3.153 

0.299 

0.345 

1.39E-02 

8.74E-04 

3.356 

1.33E-02 

3.54E-03 

3.51E-04 

14 

3.256 

3.039 

0.274 

0.333 

1.34E-02 

8.39E-04 

3.278 

1.31E-02 

3.51E-03 

3.45E-04 

16 

3.203 

2.960 

0.257 

0.325 

1.31E-02 

8.17E-04 

3.223 

1.29E-02 

3.50E-03 

3.43E-04 

18 

3.166 

2.907 

0.245 

0.319 

1.29E-02 

8.03E-04 

3.183 

1.29E-02 

3.50E-03 

3.41E-04 

20 

3.140 

2.872 

0.236 

0.315 

1.27E-02 

7.94E-04 

3.157 

1.32E-02 

3.52E-03 

3.41E-04 

22 

3.125 

2.850 

0.229 

0.313 

1.27E-02 

7.91E-04 

3.140 

1.41E-02 

3.54E-03 

3.43E-04 

24 

3.116 

2.839 

0.224 

0.311 

1.27E-02 

7.91E-04 

3.131 

1.63E-02 

3.57E-03 

3.47E-04 

26 

3.114 

2.836 

0.221 

0.310 

1.28E-02 

7.94E-04 

3.129 

2.04E-02 

3.62E-03 

3.53E-04 

30 

3.124 

2.848 

0.217 

0.310 

1.30E-02 

8.06E-04 

3.141 

3.84E-02 

3.82E-03 

3.82E-04 

40 

3.204 

2.948 

0.218 

0.313 

1.39E-02 

8.66E-04 

3.228 

1.78E-01 

5.29E-03 

6.11E-04 

50 

3.327 

3.102 

0.227 

0.319 

1.51E-02 

9.48E-04 

3.360 

4.83E-01 

8.52E-03 

1.12E-03 

60 

3.474 

3.283 

0.241 

0.325 

1.64E-02 

1.05E-03 

3.512 

9.10E-01 

1.32E-02 

1.84E-03 

80 

3.798 

3.687 

0.275 

0.335 

1.93E-02 

1.27E-03 

3.841 

1.86E-F00 

2.39E-02 

3.49E-03 

100 

4.137 

4.115 

0.315 

0.342 

2.24E-02 

1.52E-03 

4.178 

2.66E-b00 

3.34E-02 

4.95E-03 

120 

4.478 

4.549 

0.359 

0.347 

2.54E-02 

1.81E-03 

4.514 

3.21E-b00 

4.07E-02 

6.05E-03 

140 

4.815 

4.982 

0.405 

0.351 

2.85E-02 

2.11E-03 

4.844 

3.54E-b00 

4.59E-02 

6.82E-03 

160 

5.145 

5.410 

0.452 

0.353 

3.16E-02 

2.45E-03 

5.168 

3.72E-F00 

4.95E-02 

7.35E-03 

180 

5.469 

5.830 

0.500 

0.354 

3.47E-02 

2.80E-03 

5.485 

3.78E-F00 

5.20E-02 

7.72E-03 

200 

5.783 

6.241 

0.550 

0.354 

3.77E-02 

3.17E-03 

5.793 

3.77E-F00 

5.38E-02 

7.99E-03 

250 

6.528 

7.218 

0.674 

0.352 

4.50E-02 

4.18E-03 

6.525 

3.58E-F00 

5.65E-02 

8.45E-03 

300 

7.204 

8.114 

0.796 

0.349 

5.18E-02 

5.23E-03 

7.193 

3.29E-b00 

5.81E-02 

8.81E-03 
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TABLE III. Low temperature rotational de-excitation thermal rate coefficients kij^i/ji(T) in the 
HD(i) -|- ortho-ll 2 {j) —)• HD(i') -|- ortho-ll 2 {j') collision. All results are multiplied by a constant 
value a = 10^^. The data are in the unit cm^ sec“^. 


r(K) 

11-01 

21-11 

21-01 

31-21 

31-11 

31-01 

03-21 

03-11 

03-31 

13-03 

13-31 

2 

4.921 

4.903 

0.322 

3.691 

0.299 

3.68E-02 

6.98E-02 

7.74E-03 

2.16E-03 

5.037 

1.50E-01 

4 

4.189 

3.807 

0.256 

2.669 

0.218 

2.69E-02 

3.90E-02 

4.30E-03 

1.20E-03 

4.222 

l.lOE-01 

6 

3.844 

3.383 

0.229 

2.304 

0.189 

2.33E-02 

2.99E-02 

3.28E-03 

9.19E-04 

3.855 

9.58E-02 

8 

3.623 

3.146 

0.214 

2.112 

0.173 

2.14E-02 

2.57E-02 

2.80E-03 

7.94E-04 

3.625 

8.88E-02 

10 

3.472 

2.997 

0.204 

1.996 

0.164 

2.03E-02 

2.34E-02 

2.54E-03 

7.31E-04 

3.469 

8.48E-02 

12 

3.364 

2.898 

0.198 

1.922 

0.158 

1.96E-02 

2.20E-02 

2.38E-03 

7.00E-04 

3.359 

8.24E-02 

14 

3.287 

2.831 

0.194 

1.874 

0.154 

1.91E-02 

2.12E-02 

2.28E-03 

6.92E-04 

3.280 

8.10E-02 

16 

3.231 

2.787 

0.191 

1.843 

0.152 

1.88E-02 

2.07E-02 

2.22E-03 

7.02E-04 

3.223 

8.03E-02 

18 

3.192 

2.759 

0.189 

1.825 

0.151 

1.87E-02 

2.04E-02 

2.18E-03 

7.30E-04 

3.183 

8.00E-02 

20 

3.164 

2.743 

0.189 

1.816 

0.150 

1.87E-02 

2.03E-02 

2.17E-03 

7.77E-04 

3.155 

8.00E-02 

22 

3.147 

2.737 

0.189 

1.814 

0.150 

1.87E-02 

2.03E-02 

2.16E-03 

8.42E-04 

3.137 

8.04E-02 

24 

3.138 

2.737 

0.189 

1.817 

0.151 

1.88E-02 

2.04E-02 

2.17E-03 

9.26E-04 

3.127 

8.09E-02 

26 

3.134 

2.744 

0.190 

1.826 

0.152 

1.89E-02 

2.06E-02 

2.19E-03 

L03E-03 

3.124 

8.16E-02 

30 

3.143 

2.772 

0.193 

1.853 

0.155 

1.94E-02 

2.11E-02 

2.25E-03 

1.28E-03 

3.131 

8.33E-02 

40 

3.220 

2.895 

0.205 

1.962 

0.166 

2.10E-02 

2.30E-02 

2.47E-03 

2.12E-03 

3.207 

8.87E-02 

50 

3.343 

3.060 

0.220 

2.105 

0.181 

2.31E-02 

2.54E-02 

2.78E-03 

3.13E-03 

3.328 

9.50E-02 

60 

3.488 

3.249 

0.238 

2.269 

0.199 

2.56E-02 

2.80E-02 

3.14E-03 

4.19E-03 

3.471 

1.02E-01 

80 

3.811 

3.661 

0.279 

2.634 

0.240 

3.16E-02 

3.35E-02 

3.98E-03 

6.40E-03 

3.792 

1.14E-01 

100 

4.150 

4.094 

0.323 

3.028 

0.287 

3.88E-02 

3.92E-02 

4.96E-03 

8.72E-03 

4.127 

1.26E-01 

120 

4.490 

4.532 

0.370 

3.440 

0.339 

4.69E-02 

4.50E-02 

6.07E-03 

1.12E-02 

4.465 

1.37E-01 

140 

4.827 

4.969 

0.419 

3.864 

0.396 

5.60E-02 

5.08E-02 

7.28E-03 

1.38E-02 

4.799 

1.46E-01 

160 

5.157 

5.400 

0.470 

4.295 

0.457 

6.60E-02 

5.66E-02 

8.58E-03 

1.65E-02 

5.127 

1.55E-01 

180 

5.480 

5.824 

0.521 

4.731 

0.521 

7.69E-02 

6.24E-02 

9.97E-03 

1.93E-02 

5.447 

1.62E-01 

200 

5.795 

6.238 

0.573 

5.169 

0.589 

8.87E-02 

6.80E-02 

1.14E-02 

2.22E-02 

5.759 

1.69E-01 

250 

6.539 

7.224 

0.704 

6.258 

0.772 

1.21E-01 

8.14E-02 

1.53E-02 

2.98E-02 

6.498 

1.82E-01 

300 

7.215 

8.127 

0.831 

7.311 

0.964 

1.57E-01 

9.37E-02 

1.93E-02 

3.73E-02 

7.171 

1.92E-01 
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FIG. 1. (Color online) Four-atomic system (12)-|-(34) or HD-I-H 2 . Here, H is a hydrogen atom 
and D is deuterinm, represented by the few-body Jacobi coordinates: Ri, R 2 and R 3 . The vector 
connects the center of masses of the HD and H 2 molecules, i.e. O and Oi respectively, and 
is directed over the axis OZ, 61 is the angle between Ri and R 3 , 62 is the angle between R 2 and 
J? 3 , (/?2 is the torsional angle, ji, j 2 and L are quantum angular momenta over the corresponding 
Jacobi coordinates Ri, R 2 and iis. 
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FIG. 2. (Color online) The total cross section of the HD(1) + H2(0) —)■ HD(0) + H2(0) inelastic 
rotational energy transfer collision. The numbers in the brackets are the rotational quantum 
numbers of the corresponding two-atomic molecules. 
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FIG. 3. (Color online) The total cross section of the HD(2) + H2(0) —)■ HD(1) + H2(0) inelastic 
rotational energy transfer collision. The numbers in the brackets are the rotational quantum 
numbers of the corresponding two-atomic molecules. 
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FIG. 4. (Color online) The total cross section of the HD(2) + H2(0) —)■ HD(0) + H2(0) inelastic 
rotational energy transfer collision. The numbers in the brackets are the rotational quantum 
numbers of the corresponding two-atomic molecules. 
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FIG. 5. (Color online) Upper plot: the rotational de-excitation thermal rate coefficients kij^i/j' (T) 
for the HD(1) -|-H2(0) ^ HD(0) -|-H2(0) collision. Middle and lower plots represent kij^i'ji(T) for 
the HD(2)-|-H2(0) —> HD(1)-|-H2(0) and HD(2)-|-H2(0) HD(0)-|-H2(0) channels correspondingly. 
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FIG. 6. (Color online) Upper plot: the total cross section of the HD(0) + H2(2) HD(2) + 

H2(0) inelastic rotational energy transfer collision. The numbers in the brackets are the rotational 
quantum numbers of the corresponding two-atomic molecules. Lower plot: the corresponding 
thermal rate coefficient, i.e. the process (02)-(20). 
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Relative Velocity v (m / sec) 


FIG. 7. (Color online) The total cross section (upper plot) and thermal rate coefficients (lower 
plot) of the HD(1) + H2(2) —> HD(0) + H2(2) inelastic rotational energy transfer collision. The 
numbers in the brackets are the rotational quantum numbers of the corresponding two-atomic 
molecules. 
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FIG. 8. (Color online) The total cross section of the HD(1) + H2(2) —)■ HD(2) + H2(0) inelastic 
rotational energy transfer collision. The numbers in the brackets are the rotational quantum 
numbers of the corresponding two-atomic molecules. 
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FIG. 9. (Color online) The total cross section of the HD(1) + H2(2) —)■ HD(1) + H2(0) inelastic 
rotational energy transfer collision. The numbers in the brackets are the rotational quantum 
numbers of the corresponding two-atomic molecules. 
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FIG. 10. (Color online) The total cross section of the HD(1) + H2(2) —> HD(0) + H2(0) inelastic 
rotational energy transfer collision. The numbers in the brackets are the rotational quantum 
numbers of the corresponding two-atomic molecules. 








Relative Velocity v {rr\ / sec) 



Relative Velocity v (m / sec) 


FIG. 11. (Color online) Upper graph: the total cross sections of the HD(2) + H2(2) —> HD(1) + 
H2(2) and HD(2) + H2(2) ^ HD(3) + H2(0) inelastic rotational energy transfer collisions. Lower 
graph: the total cross sections of the HD(2) + H2(2) —> HD(0) + H2(2) and HD(2) + H2(2) —>■ 
HD(2) + H2(0) inelastic rotational energy transfer collisions. The numbers in the brackets are 
the rotational quantum numbers of the corresponding two-atomic molecules. Results are obtained 
with two different PESs. 
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FIG. 12. (Color online) The HD+ori/io-H 2 case. Upper plot: the total cross section of the HD(1) 
+ H2(1) ^ HD(0) + H2(1) inelastic rotational energy transfer collision. Middle and lower plots 
represent the following two integral cross sections: HD(2) + H2(1) ^ HD(1) + H2(1) and HD(2) + 
H2(1) —7> HD(0) + H2(1) correspondingly. The numbers in the brackets are the rotational quantum 
numbers of the corresponding two-atomic molecules. 
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FIG. 13. (Color online) The HD+ort/io-H 2 case. Upper plot: the rotational de-excitation thermal 
rate coefficients kij^i/j'{T) for the HD(1) -|- H 2 (l) —)■ HD(0) -|- H 2 (l) collision. The middle and 
lower plots represent kij^iiji[T) for the HD(2) -|- H 2 (l) —)• HD(1) -|- H 2 (l) and HD(2) -|- H 2 (l) —>■ 
HD(0) -|-H2(1) channels correspondingly. 
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FIG. 14. (Color online) The HD+ort/io-H 2 case. Upper plot: the rotational de-excitation ther¬ 
mal rate coefficients for the HD(0) -|- H2(3) —>■ HD(2) -|- H 2 (l) and HD(0) -|- H2(3) —>■ 

HD(1) -|- H2(1) collisions together with the corresponding result froroii. The lower plot represents 
for the HD(0) + H2(3) ^ HD(3) + H2(1) and HD(0) + H2(3) ^ HD(0) + H2(1) channels 
correspondingly. The results of this paper were obtained with the use of the modihed Hinde and 
BMKP PESs. 
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FIG. 15. (Color online) The HD+ort/io-H 2 case. Upper plot: the rotational de-excitation thermal 
rate coefficients kij^i'ji{T) for the HD(1) -|- H2(3) —)• HD(0) -|- H2(3) collision. The middle plot 
represents kij^iijiiT) for the HD(1) -|- H2(3) ^ HD(3) -|- H 2 (l) process, and lower plot shows 
three results for the HD(1) + H2(3) ^ HD(2) + H 2 (l), HD(1) + H2(3) ^ HD(1) + H 2 (l), and 
HD(1)-|-H2(3) ^ HD(0)-|-H2(1) collisions correspondingly. The results of this paper were obtained 
with the use of the modified Hinde and BMKP PESs. 
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FIG. 16. (Color online) Four-atomic system 1234 or H-D-H-H is shown together with its few-body 
Jacobi coordinates {Ri, R 2 , Rs}- The original cartesian coordinate system is OXYZ. The center of 
mass of the original H 2 molecule lies in point O. O'X'Y'Z' is a system which was shifted in parallel 
from the original system, O' lies in the center of mass of the actual HD molecule. The close-coupling 
equations are solved using the space-hxed coordinate system O'X"Y" Z". The vector R'^ connects 
the center of masses of the HD and H 2 molecules. HD is the first molecule with a rotational 
constant Be(1)=44.7 cm“^ and quantum angular momentum of ji, H 2 is the second molecule with 
rotational constant of Be(2)=60.8 cm“^ and uses momentum j 2 in the system; Ri = 0.7631 A and 
i?2 = 0.7668 A are fixed interatomic distances in each hydrogen molecule HD and H 2 respectively. 
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